The Supernova Remnant Gil. 2— 0.3 and its central Pulsar 
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ABSTRACT 

The plcrion inside the composite Supernova Remnant Gil. 2-0. 3 appears to 
be dominated by the magnetic field to an extent unprecedented among well 
known cases. We discuss its evolution as determined by a central pulsar and the 
interaction with the surrounding thermal remnant, which in turn interacts with 
the ambient medium. We find that a plausible scenario exists, where all the 
observations can be reproduced with rather typical values for the parameters of 
the system; we also obtain the most likely period for the still undetected pulsar. 

Subject headings: ISM: Supernova Remnants — ISM: Individual: Gil. 2-0. 3 — 
Stars: Pulsars: General 
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1. Introduction 

Vasisht et al. (1996, hereafter VA) have recently discussed the nature of the remnant 
Gil. 2-0. 3 (Historical Supernova 386 AD) and have shown that it is likely to be of a 
composite nature, as already suggested on the basis of radio data by Morsi & Reich (1987). 
The hard X-ray emission observed by ASCA, combined with previous Einstein data and 
observations at other wavelengths, strongly suggests (despite the fact that a pulsar has not 
yet been discovered) the presence of a plerion surrounded by a thermal shell. The data 
across the electromagnetic spectrum can then be combined with the theory of the evolution 
of pulsar-powered remnants (Pacini & Salvati 1973 (PS); Bandiera, Pacini & Salvati 1984; 
Reynold & Chevalier 1984) and allow a determination of the physical parameters of the 
system. The main outcome of the interpretation proposed by VA is that the plerion 
should be characterized by a magnetic energy much higher than the energy channelled into 
relativistic electrons, at variance with, e.g., the Crab Nebula. 

In our paper we study the evolution of the plerion coupled with that of the surrounding 
shell. We show that the magnetic field determination of VA is only marginally consistent, 
from a dynamical point of view, if the thermal and non-thermal remnants are in pressure 
equilibrium. We identify, on the other hand, a scenario that is in agreement with the 
observations, and derive the initial and present parameters of the central pulsar. In 
particular, we show that the apparent preponderance of the magnetic component in the 
pulsar output may be a natural outcome of the evolution, and does not imply a strong 
unbalance in the mechanism of production. 

Following VA we assume that Gil. 2-0. 3 is the remnant of SN 386 AD and therefore has 
an age t — 1610 yr = 5 x 10^*^ s. For its distance we use a reference value of 5 kpc and define 

— d/(5kpc). The outer thermal shell has radius Rs — 3.3d5pc, electron temperature 
T, = 0.73 keV = 8.5 x 10^ K, X-ray luminosity L, - lO^^rf^ ergs'^ (in the 0.6-3.3 keV 
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band), from which VA derive an ambient density Uo ^ 1-5 cm~^, by assuming a Sedov 
(1959) expansion in the interstellar medium. The central plerionic component has a radius 
Rp Id^ pc and its spectrum is equal to 

^- " ^ ( 2.5 xlO.^Hz )""''^ <■» 

corresponding to an X-ray luminosity Lp ~ 1 x lO^^df ergs~^, in the 0.5-4 keV band. cq. (1) 
is consistent with the observations down to 32 GHz, but overestimates the plerion emission 
at 1 GHz. Accordingly, the spectrum must have a break at a frequency i/g between 1 GHz 
and 32 GHz. 



2. The model 

The assumption of the PS model is that the central pulsar feeds both magnetic field 
and relativistic particles into the expanding plerion. Taking into account the adiabatic 
losses and the rate of production, the evolution of the field strength is given by 

Here p is the fraction of the pulsar energy loss S that goes into magnetic energy; it is 
assumed that p is constant. According to standard pulsar electrodynamics Cl oc —11", with 
n being the so-called braking index, and 

^^(l+t/r„r' 

with a — {n + l)/{n — 1) = 2 in a pure dipole field. Measured values of n range roughly 
between 1.4 for the Vela pulsar (Lyne, private communication) and 2.8 for PSR 1509-58 

(Kaspi et al. 1994). The pulsar is assumed to inject also relativistic particles, whose 
energies are distributed according to a power law with constant index 7, at a rate (1 — p)S. 
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The evolution of the electron energies is determined by both adiabatic and synchrotron 
losses. 

An immediate implication of the theory is that the emitted spectrum of the plerion 
should show a break at a frequency Uh such that particles radiating at Ub have a lifetime 
equal to the age of the nebula. Since Uf, oc B~H~'^, if the observed break is identified with 
Ub one directly derives the field strength in Gil. 2-0. 3 

which corresponds to a total field energy 

The synchrotron spectrum above Ub, up to a maximum frequency Vmax, is given by 

Prom the observations 7 = 1.6 and v^ax > 2.5 x 10-*^^ Hz; by combining eqs. (1) and (6) 
we find the present-time energy loss from the pulsar: 

^35 / '^max \ "5 „„„„-! 



£ = 2.24xlO^M ^::^ergs-\ (7) 

V2.5 X 10^8 Hzy 1-p ^ ^ 

We describe the relative importance of the magnetic and particle components in terms 
of the dimensionless quantity Q — B^Rp/6{l — p)£t. By substituting eq. (3) into eq. (2) 



dx- (8) 



and integrating, we find 

p {i+ir /■* 

1-p 1^+"- Jo (l + xY 
here we put t = t/To, and assumed a power law for the expansion, Rp oc t^. Note that Q is 
a measure of the balance at production rather than a gauge of equipartition: in fact, when 
p «i 0.5, (5 ~ 1 until t < To (and later on as well if a = 1 + r). The ratio of magnetic to 
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particle energy, instead, can be arbitrarily high if the radiative hfetime of the energetically 
important particles is sufficiently shorter than t. On the observational side, eqs. (5) and (7) 
give 

the magnitude of Q in the present case, to be compared with Q ~ 1 for the Crab plerion, is 
a quantitative measure of the discrepancy which we want to address. For Gil. 2-0. 3, then, 
either p is very close to unity (1 — p < 10~^), or f 3> 1, {a — 1 — r) > 0. We regard the 
former possibility as very unlikely, and in the following we shall limit ourselves only the 
latter. 

When t ^ To, it is appropriate to assume flux conservation for the evolution of the 
magnetic field, B oc t~^^; furthermore, the plerion spectrum shows two breaks (PS), of 
which the higher one, Ufj, corresponds to the radiative lifetime at the current time t, and 
the lower one, Vc, is the signature of the radiative lifetime at Tq. If we identify the observed 
break with Uf,, eq. (4) holds, and the pressure interior to the plerion is relatively very high. 
As we show in the following, the magnetic pressure corresponding to eq. (4) exceeds the 
shell pressure by a large factor; this is possible only if the plerion expansion is limited by 
the inertia of some cold matter tied to the plerion itself, in which case we expect r > 1. 
Then r > 1 and a < 3, as observed in all measured cases, give p of order 0.5 only at the 
expense oi t > 1000, and an initial magnetic energy of order i x Wb > lO^^dg erg. We 
consider very unlikely a value so much larger than the canonical supernova output, and 
propose that: i) the plerion and the shell are in pressure equilibrium, so that r < 1; and 
ii) the observed break is to be identified with i/c, so that — Uc x ^ i/g, and Q 

is decreased with respect to eq. (9). However, the break at uj, entails a change of slope 
e = (2.8r — 0.4Q;)/(5r — 1) (see PS); since Gil. 2-0. 3 shows a straight spectrum from the 
microwaves to the X rays, we must require e = 0, a = 7r. 
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Within this framework more detailed results can be obtained by coupling the plerion 
evolution to the hydrodynamics of the thermal shell. We can distinguish two different cases, 
according to the distribution of the ambient density being of the interstellar type (po ~ 
constant, henceforth IS), or circumstellar type (po ~ henceforth CS). If the ambient 
medium has a constant density, Rg oc t^/^, Rp oc and, because of the condition on e, 
a = 2.1. If the shell expands inside the pre-supernova wind, one has Rg oc t^/^, Rp oc t^/^, 
and a — 4.0. 

We have computed in detail the two cases, and give the results in Table 1. In particular, 
we have adjusted the unknown distance 0^5 so as to make the total energy of the shell 
remnant equal to the canonical 10^^ erg. We have chosen the value of Uo within the allowed 
interval so as to bring our results as close as possible to 'typical' values. And, finally, we 
have neglected the work done by the plerion on the shell remnant, which is justified a 
posteriori by the relatively small plerion energetics resulting from the computations. The 
precise value of the ambient density in the IS and CS cases was computed so as to reproduce 
the temperature and flux of the thermal X rays given in VA. 

Note that the CS assumption leads to very 'palatable' estimates for all the unknowns, 
with one exception: the slowing down exponent a implies a braking index n = 1.67. We 
note that such value is still in the acceptable range between a purely dipolar {a — 3, n — 3) 
and a purely monopolar geometry {a — 00, n — 1), but that for most pulsars the braking 
index is between 2 and 3. On the other hand, the IS hypothesis carries along some deviation 
from perfect balance at production (p ~ 0.9) and a pulsar magnetic field at the extreme of 
the observed range {B^ = 1.3 10^^ G, like in PSR 1509-58). All in all, we regard the latter 
interpretation as the most plausible one, also because it places Gil. 2-0. 3 exactly on the 
empirical relation between Lp and S given by Seward & Wang (1988); wc suggest that the 
corresponding values in Table 1 be taken as a reasonable guess for further observing efforts. 
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A major improvement in the modelling could ensue if future observations were to give 
the precise frequency of the break(s) in the spectrum and the related change(s) of slope, or 
the detection of the pulsar. 

3. Conclusions 

If Gil. 2-03 is a composite remnant deriving from the historical explosion of SN 386 
AD, the main parameters of the nebula and of the central pulsar can be estimated in 
the framework of the PS model. They appear to be consistent with the properties of the 
remnant and the standard views on pulsars. 

The apparent preponderance of the magnetic channel in the output of the central pulsar 
is the consequence of the dynamical evolution under the influence of the external shell. 
This has slowed down the expansion and therefore reduced the importance of the adiabatic 
losses with respect to the synchrotron losses. Since most of the energy was released by the 
pulsar during a relatively short initial phase, it is natural that the subsequent evolution 
has led to a large value of Q, eq. (8), even if the pulsar output was always well balanced, 
p 0.5. We strongly stress this point as relevant to all sorts of non-thermal sources: only 
if the generation of particles and fields keeps occurring for a time scale longer or roughly 
equal to the age of the source does Q reflect the balance at production between fields and 
particles. This is why Q 1 in the Crab plerion, for which t 3 x Tq. In addition, and 
for quite independent reasons, the Crab is not far from equipartition, since its V}, is not far 
from the frequency where most of the power is emitted. 

This work was partly supported by a Grant of the Italian Space Agency. 
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Table 1. Parameters of the Models 



Parameter IS Model CS Model 



Vo 4 GHz 32 GHz 

d 9.2 kpc 7.4 kpc 

p 0.9 0.3 

Bp 9.3xlO-^G 1.8xlO-^G 

P 0.34 s 2.73 s 

£ 1.55 X lO^^ erg s"^ 1.25 x 10^6 erg s"^ 

2.5xlO^^G 1.6 X 10^2 Q 

To 5.8xl0«s 7.0xl0'^s 

Po 0.029 s 0.139 s 

£o 2.4 X lO^'^erg 1.0 x 10^^ erg 
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